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Abstract 

A two-dimensional Schrodinger operator with a constant magnetic field perturbed by a 
smooth compactly supported potential is considered. The spectrum of this operator consists 
of eigenvalues which accumulate to the Landau levels. We call the set of eigenvalues near the 
n'th Landau level an n'th eigenvalue cluster, and study the distribution of eigenvalues in the 
n'th cluster as n — > oo. A complete asymptotic expansion for the eigenvalue moments in the 
n'th cluster is obtained and some coefficients of this expansion are computed. A trace formula 
involving the first eigenvalue moments is obtained. 



1 Introduction and Main Results 

1. Introduction Let H in L 2 (R 2 , dx\ dx<i) be the following magnetic Schrodinger operator: 

d B \ 2 ( . d B 2 



The operator H describes a quantum particle in R 2 in a constant homogeneous magnetic field of 
the magnitude B; it is often called the Landau Hamiltonian. It is well known [8] that the spectrum 
of H consists of a sequence of eigenvalues (Landau levels) A n = B(2n + 1), n E Z + = {0, 1, 2, . . . }. 
Each of these eigenvalues has infinite multiplicity. 

Let V G C^°(R 2 ) be a real valued function (potential in physical terminology). Consider 
the spectrum of the operator H + V. It is well known (see [3]) that V is a relatively compact 
perturbation of H and therefore the essential spectrum of H + V is the same as that of H, i.e. 
consists of the Landau levels. The operator H + V may have eigenvalues of finite multiplicities 
which can accumulate to the Landau levels. 

Let us define disjoint intervals Ao = [inf a(H+V), 2B), A n = [A n — B, A n + B), n G N, so that 
a(H + V) C U^qA^. We shall call the set a(H + V) n A n the n'th eigenvalue cluster. For a fixed 
n, the distribution of eigenvalues in the n'th cluster was studied in [10], [11] (these papers contain 
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also references to earlier work on this problem). It was found that eigenvalues accummulate to 
A n exponentially fast. 

Our aim is to study the asymptotic distribution of eigenvalues in the n'th cluster as in oo. 
Our first preliminary result is that the width of the n'th cluster is 0(n~ 1 / 2 ): 

Proposition 1.1. There exist C > and N G N such that for all N, one has 

a(H + V) n A n C (A n - Cn- 1 / 2 , A n + Cn" 1 / 2 ). 

The constants C and N depend only on sup|^(x)| and on the diameter of suppV . 

xem. 2 

The power n -1 / 2 in the above Proposition is sharp; see Remark 3.2 below. 

2. Definition of eigenvalue moments fi n . We would like to define moments of eigenvalues 
in the n'th cluster. First, in order to explain the main idea of the definition, let us define the 
eigenvalue moments 'naively' for the case ||V|| < B; here and in what follows ||V|| = ||V||loo. 
Fix n G Z + and enumerate Ai, A2, A3,. . . all eigenvalues in the n'th cluster so that |Ai — A n | ^ 
I A2 — A n I ^ I A3 — A„ I ^ • • • . Let us define the eigenvalue moments by 

= X>,' - A n ), n G 1+ (\\V\\ < B). (1.1) 
j 

By the above quoted result of [10], [11], the rate of convergence Aj — ► A n as j — ► 00 is exponential, 
and therefore the series (1.1) converges absolutely. 

In order to give the definition of eigenvalue moments which is suitable both for the case 
||V|| < B and for the case ||V|| ^ B, we need to recall the notion of the spectral shift function 
for the pair of operators H + V, H. The spectral shift function was introduced in an abstract 
operator theoretic setting in [9, 7]; see also the book [14]. Recall that (see [3]) 

(H + V - Ao)" 1 - (H - Ao)- 1 G Trace class, A < inf a(H + V). (1.2) 

This enables one to define the spectral shift function £ G L 1 1 oc (IR) for the pair H + V, H. The 
spectral shift function £ is determined by the following two conditions: 

(i) For any 'test function' (f> G Cg°(]R), one has the trace formula: 

/oo 
£W(A)dA. (1.3) 
-00 

(ii) £(A) = for A < inf a{H + V). 

Note that 4>{H + V) — 4>(H) G Trace class by (1.2) (see [14]). In fact, the class of admissible test 
functions <p is much wider than Cg°(R). In particular, this class includes exponentials 4>{\) = e _ ' A , 
t > 0; we will use the latter fact in the sequel. 

Condition (i) determines the spectral shift function up to an additive constant; condition (ii) 
fixes this constant. As it follows from the trace formula (1.3), for A G M \ a(H + V) the spectral 
shift function can be determined by (see [14], section 8.7 and formula (8.2.20)) 

£(\) = Tr{E H {\)-E H+v {\)), \€R\a(H + V), (1.4) 

where Eh(X) and Eh+vW are the spectral projections of H and H + V associated with the 
interval (—00, A). In particular, it follows that £ is constant and integer- valued on the intervals of 
the set R\a(H + V). 
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Now we are ready to give a general definition of the eigenvalue moments: 

Hn = I C(A)dA, neZ+. (1.5) 

Let us explain why the definitions (1.5) and (1.1) coincide for ||V|| < B. From (1.4) one can see 
that for II VII < B 



the number of eigenvalues of H + V in (A, A n + B) if A G (A n , A n + B); 
(—1) x the number of eigenvalues of H + V in (A n — B, A) if A G (A n — £?, A n ). 

(1.6) 

From here it follows that (1.5) and (1.1) coincide. 
Remark. Proposition 1.1 shows that 

a(H + tV) n A n C (A n - Cn" 1 / 2 , A n + Cn~ 1/2 ), Vr G [0, 1] 

for any V and all sufficiently large n. From here, using (1.4) and a continuity in t argument, one 
can prove that for any V and all sufficiently large n, 

supp£ n A n C (A n - Cn- 1 ' 2 , A n + Cn^ 2 ) (1.7) 

and (1.6) holds true. Thus, definition (1.1) is applicable for any V and all sufficiently large n. 
3. Main result 

Theorem 1.2. The asymptotic expansion 



"1 "3 /-. n\ 

/U n ~Q + — r-^ H ^^T^H , ra^oo, (1.8) 



holds true with some real coefficients ctj . Moreover, one has 

a o = 7T [ V(x)dx, ai = a 2 = 0, a 3 = - —^t [ [ Y^XM ' dxdy. (1.9) 
2vr J R 2 16V2tt 3 Jr2 J K 2 \x-y\ 

The identity (trace formula) 
holds true. 

Remarks (1) The coefficients ceo, «2 are obtained by comparing the asymptotic expansion 
(1.8) with the small t asymptotic expansion of Th:(e~ tlyH+v ^ — e~ tH ) — see section 2 below. It 
does not seem possible to obtain the coefficient 03 by using a similar argument; we obtain it 
by a more direct analysis (see end of section 4). (2) Some results concerning the eigenvalue 
distribution in clusters for large n can be found in [12]. (3) Similar trace formula for the two- 
dimensional perturbed harmonic oscillator was obtained in [6]. (4) It might be interesting to note 
in connection with the formula for 03 that the integral j J —^^^-dxdy appears as the coefficient 
of the leading term in the high energy asymptotic expansion of the total scattering cross-section 
for the pair of operators —A, —A + V(x) in L 2 ' 
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Along with the moments fi n , we will use the higher order moments 

M <*) = (fc + i) f (A — A n ) k £(\)d\, ken, neZ + . (1.11) 

In the case ||F|| < B, the last definition becomes ^ = ~ ^n) k+1 - We will also prove the 

asymptotic expansion 



^ k) -n'^(a^ + % + ^ + % + ---), fcGN, n ^ oo. (1.12) 

Below for consistency we write [i n = [iffl , ctj = a~p . 

4. The structure of the paper We will use three distinct arguments to prove Propo- 
sition 1.1, the asymptotic expansions (1.8), (1.12) and the trace formula (1.10). The proof of 
Proposition 1.1 is self-contained, the proof of the asymptotic expansions (1.8), (1.12) depends 
on the estimates (3.2), (3.3) obtained in the proof of Proposition 1.1, and the proof of the trace 
formula (1.10) depends on both Proposition 1.1 and the expansions (1.8), (1.12). 

In section 2, assuming Proposition 1.1 and the existence of the asymptotic expansions (1.8), 
(1.12), we prove the trace formula (1.10) and derive formulae (1.9) for the coefficients ceo, ct\ and 
a.2. The argument is quite elementary. 

Proposition 1.1 is proven in section 3. 

In Sections 4, 5, 6 we justify the asymptotic expansions (1.8), (1.12). The proof is based on a 
detailed analysis of the properties of the integral kernel of the resolvent of H (see (5.1)) and on 
some facts from the theory of confluent hypergeometric functions. This part of the paper is fairly 
elementary in nature but technically is rather complicated. 

In Section 4 we also prove the formula (1.9) for the coefficient 03. 

5. Notation We use notation \\A\\, ||^4||s 2 , \\A\\s 1 for the operator norm, the Hilbert-Schmidt 
norm, and the trace class norm of an operator A. By C, c we denote various constants in the 
estimates. 



2 Proof of the trace formula 

1. Heat kernel asymptotics. 
Lemma 2.1. The asymptotic formula 

Ti(e- tH - e ^ H+v) ) = / V{x)dx - — f V 2 (x)dx + 0{t 2 ), t -► +0 (2.1) 

4vr Jr2 8tt J R 2 

holds true. 

Proof. The required asymptotics can be obtained by using general results on asymptotic expan- 
sions of heat kernels of second order elliptic operators. However, the two term asymptotic formula 
(2.1) is considerably simpler than the aforementioned general results, and so we prefer to give a 
direct 'elementary' proof. We use the formula 

e -tH _ e -t(H+V) = t e - { t- tl )H Ve -t l{ H + V) dti (2 2) 

JO 
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and the explicit formula for the integral kernel of e tH (see [3]): 

e~ tH (x, y) = 47rgii ^ (£?t) exp(-f \x - y\ 2 coth(Bt) + i§ [x, y]), x,y £ M. 2 , t > 0, (2.3) 
where [x,y] = x\y2 — x 2 y±. Iterating (2.2), we obtain 

Tr{e -tH _ e -t(H+V)) = /l(t ) + h{ t) + h ( t ), 

I^t) = f Tv(e- ( - t ' t ^ H Ve- tlH )dt 1 = tTi(Ve- tH ) = [ V(x)dx + 0{t 2 ), t +0, 

JO 47T J R 2 

I 2 (t)= [ dh f 1 dt 2 Tr{e-^- t ^ H Ve-^- t ^ H Ve- t2H ), 
Jo Jo 

I 3 (t)= f'dh t dt 2 r^ 3 'IV(e-( t - tl ) ff F e -( tl - t2 ) ff y e -( t2 - t3 ) if Fe- t3 ( if+y )). 
Jo Jo Jo 

Consider the term I 2 (t). Introducing the new variable s = t — t\ + t 2 , we have 

I 2 (t)= f dssTx{e- sH Ve- {t - s ^ H V) 
Jo 

= / ds — — ^— — -/ dz exp(-f (coth5s + coth5(t-s))|z| 2 ) W(z), 

J smhBssmhB(t- s) J R2 FV 4V v ,n 1 ; v h 

where W{z) = / R2 + z)dy. Let us write W(z) = W(0) + (W(z) - W(0)) and split 

/2(t) accordingly: I 2 (t) = I 2 \t) + I^\t). For I^Wi explicitly computing the integrals, we 
obtain 

(2) 

For i£ ; (t), using the estimate \W(z) — W(0)| ^ C|z|, we get 

|4 2) (i)l <C /" — — ^— / dz |z|exp(-f (coth 5s + coth £(i - s))|z| 2 ) 

Jo smh B s smh B (t — s) J 4 



Ci 



sinh Bt 



f ds s (smh B s smh B(t- s)) 1/2 = 0{t 2 ) t -> +0. 
Jo 



Finally, let us estimate the term Is(t). Using the Hilbert-Schmidt norm estimate ||e tjf/ || ^ ^ 
Ct~ 1/2 , t > 0, we obtain 



|Tr(e _( * _ * l)if ye _(tl_t2)if ye _(t2_t3)H Fe _ * 3(jff+y) 



\/t — — t 2 

which yields J 3 (t) = 0(t 2 ), t -> +0. ■ 

2. Auxiliary estimate 
Lemma 2.2. One /ias 

/ |£(A)|dA = O(l), 



n — > oo. 
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Proof. Recall that the spectral shift function is monotone with respect to the perturbation V. 
I.e., denoting temporarily by £(A; V) the spectral shift function corresponding to the potential V, 
we have 

if Vi < V 2 , then f (A; Vi) < f(A; V 2 ) a.e. A G R. 
Also, we have f (A; — V) = -f (A; V). Let us choose Vi, V 2 G C£°(IR 2 ) such that 

Vi > 0, V 2 > and - V 2 < V < Vi. 

Then 

-e(A;V 2 )=^(A;-V 2 )^e(A;V)^e(A;Vi), e(A;Vi)>0, £(A;V 2 )^0. 
Therefore, |£(A; V)| < £(A; Vi) + £(A; V 2 ) and 

/ \i(\-V)\d\^ f i(\-V 1 )d\+ f £(A;V 2 )dA. 
ja„ ja„ ja„ 

By the asymptotics (1.12) for fin\ applied to the potentials V\ and V 2 , the r.h.s. of the last 
inequality is O(l) as n ^ oo. ■ 

3. Proof of the trace formula (1.10) and formulae (1.9) for • 
1. By Krein's trace formula (1.3) with 0(A) = e~ tx , t > 0, and Lemma 2.1, we obtain 

£(\)e- tx d\ = - Ti(e~ tH - e - t(H+v y) = -i- f V{x)dx - — [ V 2 {x)dx + 0(t), (2.4) 
J-oo t 4trt J R 2 8tt J R 2 

as t — > +0. On the other hand, one can rewrite the integral in the l.h.s. of (2.4) as a sum over 
the eigenvalue clusters: 

/oo OO „ 

£(A)e-* A dA = J2 £(A)e-' A dA. (2.5) 

Let us use Taylor's formula for e" tA , A G A n : 

|e- u -e-* A «(l-i(A-A n ) + ^ 2 (A-A n ) 2 )KCt 3 |A-A n | 3 , A G A n . (2.6) 
By (1.7) and Lemma 2.2, we have 

oo „ oo 

W |£(A)||A-A n | 3 dA = ^0(n- 3 / 2 )<oo. (2.7) 

n=0 An n=0 

Combining (2.5)— (2.7), we obtain 

/qq OO OO ^ oo 

«A) e -*\iA = £ ^e-* A « - i £ ^e"*^ + > £ ^-tA. + 0(t s } > t ^ +Q _ (2 . 8) 

-°° n=0 n=0 n=0 

Below we compare (2.4) and (2.8). 
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2. We will use the following elementary formulae for t — > +0: 

e- m 1 
1 - e~ 2B * ~ 2~Bt 



00 -Bt 



- p—tst i 

n=0 



/-oo /— 

E ^" 1/2 e-* A " = _/ x-^e-^^+^dx + 0(1) = + o(l); (2.10) 

00 /'OO 

Y j n- 1 e~ thn = x- 1 e- tB{2x+1) dx + 0{l) = -\ogt + 0{l); (2.11) 
n=i 71 

OO pi CO 

Vn- 3 / 2 (l-e-* A ") = / (Vn- 3 / 2 A n e- sA ")ds = 2\/2^Bt + 0(t). (2.12) 

n=l J ° n=l 

Using (2.9)-(2.11) and the asymptotic expansions (1.8), (1.12) for //^, we obtain for t — > +0: 

oo (0) (0) ^- 

E = ^ + ^ - 4 0) log* + O(i); (2.13) 



n=0 

oo 



^^) e -^=a«-^= + 0(l gt); (2.14) 



n=0 

oo 



^ M g.) e -tA n=0(]ogt) _ 



n=0 

Substituting this into (2.8) and comparing with (2.4), we find: 



a (°) = JL / vCxJefa, «S 0) = a< 0) - 0. 



Jr 2 

Thus, [in^ = offi + a^n -3 / 2 + 0(n~ 2 ), and so we get: 



oo 



X; mSV*^ = 4 0) E e " iA " + E(^ 0) - 4 0) ) + E(^ 0) - «o 0) )( e - tA " - 1) 

71=0 71=0 71=0 71=1 

(0) oo oo oo 

= ^+E(/# ) -4 0) )+"? ) E"-^ 

n=0 71=1 7i=l 

(0) oo 

= ffe + E(^ 0) " «o 0) ) - afhV^Bi + o(Vt), (2.15) 

n=0 

as t — ► +0. Upon comparing with (2.4), we find the trace formula (1.10) and also the formula 

a £ ) = -4Baf ) . (2.16) 
We will use (2.16) later in section 4 in order to determine the coefficient . ■ 



3 Proof of Proposition 1.1 

Let P n , n ^ 0, be the orthogonal projection onto the eigenspace of the operator H corresponding 
to the Landau level A n . An explicit formula for the integral kernel of P n is available (see e.g. 
[11]): 

P n (x,y) = ^-L n (^\x-y\ 2 )exp(-^\x-y\ 2 + i^[x,y]), x,y GR 2 , (3.1) 
where L n is the Laguerre polynomial and [x,y] = x\yi — xiy\- 
Lemma 3.1. Let V be any bounded function on M 2 with compact support. Then 

\\\V\ 1/2 P n \V\ 1/2 \\ = 0(n~^ 2 ), n^oo, (3.2) 
|||^| 1/2 ^|^| 1/2 ||5 2 = 0(^ 1/4 ), n ^ oo. (3.3) 

In the proof of Proposition 1.1, we will only need the operator norm estimate (3.2). The 
Hilbert-Schmidt norm estimate (3.3) will be used in section 4. 

Proof. We will use the following asymptotic formula for the Laguerre polynomials [5, 10.15(2)]: 

L n (t) = e*/ 2 J ( V(4n + 2)t) + R n (t), R n (t) = 0(n~ 3 / 4 ), n oo, (3.4) 
where the bound 0(n~ 3//4 ) is uniform on any bounded sub-interval of [0, oo). Let us write 

\V\ l ' 2 P n \V\ l / 2 =A n + B n , 
where A n and B n are the operators in L 2 (]R 2 ) with the integral kernels 



A n (x,y) = ^J (^T n \ x -y\y^\V{x)\ l / 2 \V{y)\ 1 / 2 , (3.5) 
B n (x,y) = |^„(f ^-yfy^^Vix^Viy)?' 2 . 
As V is bounded and compactly supported, (3.4) gives 

\\Bn\\ < ||B n || Sa = 0(n" 3 / 4 ), n - oo. (3.6) 
Let us prove the estimate (3.3). By (3.6), we only need to check that 

\\A n \\ S2 = 0(n- 1 ^). 

The latter estimate immediately follows from the explicit form of the kernel of A n and from the 
simple inequality |Jo(i)l ^ C/y/i for t > 0. 

Next, let us prove the estimate (3.2). Let A n be the operator in L 2 (M. 2 ) with the integral 
kernel 

Mx,y) = ^J (y/K\x - ylWWnViy)] 1 ' 2 . 

Note that up to a constant, A n coincides with the imaginary part of the sandwiched resolvent of 
the operator —A in L 2 (M. 2 ): 

An = — Im (\V\ 1/2 (-& " An " iOyW 2 ). 

7T 
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It is well known (see [2]) that 

Hl^l^-A - A - iOy'lV^W = 0(A -1 / 2 ), A - oo. 

Thus, we obtain 

\\A n \\=0(n- 1/2 ), n^oo. (3.7) 

Next, observe that the kernel of A n differs from that of A n by a factor e*"^'^. We are going to 
use this observation and apply the theory of 'multipliers of kernels of integral operators' [4]. Let Q 
be a sufficiently large ball in M? so that suppV C £1 and let p G L°°(Q x Q). For a Hilbert-Schmidt 
class operator T on L 2 (Q) with the integral kernel T(-, •) £ L 2 (Q x f2), let T be the operator with 
the integral kernel T(x,y)p(x,y). Evidently, T is also a Hilbert-Schmidt class operator and one 
has the estimate \\T\\s 2 ^ ||p||l°° \\T\\s 2 ■ 

Next, suppose that the mapping T ^ T sends the trace class S\ into itself and there is a trace 
class norm bound II^H^ ^ C^HTH^. Then, by duality between the trace class S\ and the class 
B(L 2 (S7)) of all bounded operators on L 2 (Q), the mapping T ^ T can be extended onto B(L 2 (fi)) 
and the norm bound ||T|| ^ C(/9) ||T|| holds true. In this case p is called a bounded multiplier on 
the class B(L 2 (ft)). 

A sufficient condition (see [4]) for p to be a bounded multiplier on B(L 2 (f2)) is 

sup||p(av)||ffs(Q) < oo, s > 1, 
x<=n 

where H S (Q) is the standard Sobolev class. Clearly, p(x, y) = e l "2 ^ satisfies the above condition, 
and therefore 

\\A n \\^C\\A n \\=0{n- 1 l 2 ), rwoo, 

which yields (3.2). ■ 

Proof of Proposition 1.1. The proof is valid for any bounded compactly supported potential. 
By the Birman-Schwinger principle, it suffices to show that for some C > and all sufficiently 
large n, 

|||Vf /2 i?(A)|F| 1/2 || < I, forallAeA n , |A-A n |>-^=, (3.8) 



>n 

where R(X) = (H - A) -1 . Choose / G N sufficiently large so that ||F||/Az < 1/2, and write R(X) 
as 

n+l „ 



u 1 Afc ~~ A 

Then, for A G A n , 

\\\v\^ 2 R(x)\vn\^ g ll|y| ^ l ^ |1/2|l +|||y| 1 ^(A)|y|^||, 

k=n—l ^ 

By the choice of I, one has IIIFI 1 / 2 ^^)!^! 1 / 2 !! < 1/2. On the other hand, by Lemma 3.1, 

i^Pkivn ^ (a + 1)Q(w -i /2) max 

\A k — A| n-l^k^n+l 



k=n— I 



Thus, we get (3.8) for sufficiently large C > 0. 
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Remark 3.2. 1. The operator norm estimate (3.2) is sharp, i.e. for any V not identically zero, 
one can prove that 

w\v\y 2 p n \vn^c/v^ 

for some c > and all sufficiently large n. Indeed, without the loss of generality assume that 
supp V contains an open neighbourhood of zero and let V be a spherically symmetric potential, 
V(x) = v(\x\) < \V(x)\. Then by Lemma 3.3 of [11], the eigenvalues A fc of \V\ l / 2 P n \V\ 1 / 2 are 
given by 

+ J v(^/2t/B)e~H k L^\t) 2 dt, k = -n,-n + l,-n + 2,..., (3.9) 

where are the Laguerre polynomials. Taking k = and using the asymptotics (3.4) and the 
asymptotics of the Bessel function, one obtains that (3.9) has asymptotic behaviour cn _1 / 2 (l + 
o(l)), as n — > oo. 

2. Using the above observation and an argument similar to the proof of Proposition 1.1, one 
can prove that Proposition 1.1 is sharp in the following sense. Suppose that the potential V is not 
identically zero and is either non-negative or non-positive. Then for some c > and all sufficiently 
large n, (the width of the n'th eigenvalue cluster) ^ cn -1 / 2 . 

4 Proof of the asymptotic expansions (1.8), (1.12) 

We will prove the asymptotic expansion (1.12) by expressing the eigenvalue moments as 
contour integrals of an analytic function. Let T n be a positively oriented circle around A ra with 
the radius B. First, we need estimates on the norm of the 'sandwiched resolvent' of H on the 
contours Y n . 

Lemma 4.1. For n — > oo ; one has 

sup|||V| 1/2 i2(z)|V| 1/2 || =0(n- 1 / 2 logn), (4.1) 

suplHVf/^O^lVf/ 2 !^ = 0(n -1 / 4 log n). (4.2) 

zer„ 

Proof. Let us prove (4.1). Using the estimate (3.2), we get for z € F n : 

00 IIIVIl/2 P |T/|1/2|| 00 r 

^ \ A k ~ A| ^ v^Afc - *| 

r l —— ^— — r +c r ^ , dx , — r +o(n- i /2 )=0(n -i/2 logn) 

7 ^\B(2x + 1) - z| 7 n+1 v ^| J B(2x + 1) - z| V ; V ; ' 
as n — ► oo. The estimate (4.2) can be proven in a similar fashion by using (3.3). ■ 

The core of the proof of the expansions (1.8), (1.12) is the following Lemma. 
Lemma 4.2. For all k € Z + and all j ^ 2, the integrals 

[ Ti(VR(z)y(z-A n ) k dz (4.3) 
have an asymptotic expansion in integer powers of n~ x l 2 as n — > oo. 
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The proof of Lemma 4.2 is given in sections 5, 6. 

Proof of the asymptotic expansions (1.8), (1.12). First of all, note that it suffices to prove 
(1.8), (1.12) with some complex coefficients df^; indeed, as /j,^ are real, a posteriori the coef- 

(k) 

flcients a - are easily seen to be real. Thus, in what follows we will work with expansions with 
complex coefficients. 

By [3, Theorem 2.11], the difference of the resolvents of H + V and H belongs to the trace 
class. This enables us to define the analytic function 

W(z) = Tt{(H + V - z)- 1 -(H- z)- 1 ) = - r ,^ X \ 2 d\ z e C \ (a(H) U a(H + V)). 

J-oo (A - ZY 

The second equality in the above formula is due to Krein's trace formula (1.3). Let n be sufficiently 
large so that 

B B 

SU P p£ n A„ C [An - — , A n + — ] 

(see (1.7)). Let, as above, F n be a positively oriented circle around A n with the radius B. 
Integrating W(z)(z — A n ) k+1 over z around T n , we obtain 

- h L w ^ - ^ - L ^ h L 

= (k + l)f e(A)(A-A„) fc dA = / u(f). (4.4) 

Firstly we prove the expansion (1.12) (i.e. assume k ^ 1). Expanding the resolvent (H + V — z)~ l 
yields 

oo 

w( Z ) = Y J (-iynR(z)(vR(z)n (4.5) 

3=1 

Lemma 4.1 ensures that the series in (4.5) converges absolutely for z G T n and large n. Substituting 
the expansion (4.5) into (4.4) and subsequently integrating by parts in each term of the series, we 
obtain 

M W = (fc + l) £ M_ / Tr(VR(z)y(z-A n ) k dz, kN. (4.6) 
j=k+l ^ m Tn 

Here the summation starts from j = k + 1, as for j ^ k the integrand is analytic at z = A n 
and therefore the integral vanishes. Using Lemma 4.1, we obtain the following estimates for the 
integrals in the r.h.s. of (4.6): 

j Tr(VR{z)) j (z-A n ) k dz < B k j \\\V\ 1/2 R{z)\V\ ll Ys 2 \\\y\ 1,2 R{ z )\y\ l,2 \\ j ' 2dz 

, , , Clogn^ 2 (C\ogn\^ 2 



n i/4 j y n i/ 2 

This ensures that the series (4.6) converges absolutely (for sufficiently large n) and gives a bound 
for the remainder: 



j=N ^" 



o 



((log n^n-M/ 2 ) , n^oo. (4.7) 
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Combining Lemma 4.2 with the estimate (4.7), we obtain that the moments Hn \ k ^ 1 have an 
asymptotic expansion in integer powers of n -1 / 2 as n — ► oo. We also need to prove that first 
several terms of the expansion vanish, so that the expansion starts from the term Cn~ k t 2 . This 
can be seen as follows. Note that for j = k + 1 we can compute the integral in the series (4.6), 
which gives 



(4.8) 



a^=Tr(VP n ) k+1 + (k + l) V { — L — / Tr(VR(z)) j (z-A n ) k dz, k G Z+. 
Lemma 3.1 gives 

Tr(FP n ) fc+1 < lllFr^Pnl^l 172 !!!,!!^! 172 ^!^ 172 !!^ 1 = 0{n- k ' 2 ), n ^ oo. 

Combining this with the estimate (4.7) with TV = k + 2, we obtain /x!^ = 0{n~ k l 2 ) as n — > oo. 

Secondly we prove the expansion (1.8), i.e., the case k = 0. Here the only difference is that 
the first term in the series (4.6) is not well defined, as VR(z) is not of the trace class. However, 
this term can be written as (cf. (4.8)) Ti{\V\ l / 2 P n \V\ 1 / 2 signy), and by the explicit form (3.1) of 
the integral kernel of P n , the last expression equals ^ J V(x)dx. The rest of the argument is the 
same as for k ^ 1. ■ 

Proof of the formula (1.9) for af\ We will prove 

( i)_ B 3 / 2 f f V(x)V(y) 



a\ > = _ / / , v ' K "' dxdy. (4.9) 
4V2ir3j R 2j R 2 \x-y\ V ' 



By (2.16), this will also imply the formula (1.9) for Oi^\ 
Due to (4.8) and (4.7), we have 

= Ti(VP n ) 2 + 0((\ognfn~ l ), n - oo, 

so it suffices to prove that 



„(i) 



Tr(VP n ) 2 = %+o(n- 1 / 2 ), n ^ oo, 



with «q given by (4.9). By formula (3.1) for the integral kernel of P n , we have 
Tr(VP n ) 2 =(^-) [ V (x)V(y)L n (§\x-y\ 2 ) 2 ex V (-§\x-y\ 2 )dxdy 

= 2^2 J Ln ^ e " t h ^ tdt ' 

where h G Cq (R) is given by 



h(t)= du I dyV(y)V(y + Jjjtu), teR. 
Js 1 Jr 2 v 

By (3.4), 

/ L n (i 2 )V* h(t)tdt= J {WAn + 2) 2 h(t)tdt + 0(n~ 3 / 4 ), 



n — > oo. 
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Next, using the asymptotics of the Bessel function, we obtain 

\Mx)- v^ cos (^-f)| ^Cx-^ 2 {i + x)-\ x>0, 

and therefore 

\Jo(x) 2 - ^(cos(x- f)) 2 | ^ Cx-\l + x)-\ x>0. (4.10) 

One has 

roc poo 

^^ 2 (^s(tV^T2-l)) 2 h(t)tdt = ^ r l h(t)dt 

J " 

/•OO /"CO 

+ wk+2 cos(2tV4n + 2 - f )/i(t)<ft = ^ / /i(0^ + °(™~ 1/2 ), n - oo; 
</ " 

(4.11) 

poo 

Combining (4.10), (4.11), (4.12), and computing the integral h{t)dt yields formula (4.9). ■ 

5 Analytic properties of the resolvent R(z) 

In this section, we discuss analytic properties of the integral kernel of the resolvent R(z) = 
(H — z)" 1 and reduce the proof of Lemma 4.2 to Lemma 6.1. Our analysis is based on the 
following explicit formula for this kernel: 

Lemma 5.1. For any z G C \ cr(H), the integral kernel of the resolvent R(z) of the magnetic 
Hamiltonian H can be expressed in terms of T -function and confluent hyper geometric function 
U (a, b; C) as follows: 

R(z)(x, y) = i-r (i - U (i - 2%, 1; f|x - y\ 2 ) exp (-f \x - y\ 2 + if [x, y}) , (5.1) 
where x,y £ M 2 , x / y. 

Proof. Let us employ the integral representation [1, (13.2.5)] for the confluent hyper geometric 
function 

roo 

T(a)U(a,l;C) = e _CT r a_1 (l + r)~ a dr, 0<Reo<l, C>0 (5.2) 
Jo 

and the explicit formula (2.3) for the heat kernel [3] of the magnetic Hamiltonian H. Substituting 
(2.3) into the formula 

roo 

R(z)= / e- tH e tz dt, 
Jo 

denoting £ = f \x — y\ 2 , making the change of variable r = (coth(Bt) — l)/2 in the integral, 
and taking into account (5.2), one obtains (5.1) for — B < Rcz < B. Analytic continuation in z 
completes the argument. ■ 
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For the reader's convenience and ease of further reference we start by recalling the necessary 
facts about the confluent hyper geometric functions U(a,b;Q, M{a,b\Q. Our main sources are 
[1] and [13]. The functions U(a,b;() and M(a,b;() are two linearly independent solutions to the 
Kummer's equation 

d 2 U „ ^dU 

We are only interested in the case 6 = 1 (see (5.1)) or b lying in a small neighbourhood of 1, so we 
assume that |6 — 1| < 1/2; this will simplify our discussion. We also assume < £ ^ R for some 
fixed R > 0, as we are interested in the case £ = ^\x — y\ 2 when both x and y are in suppF (see 
(4.3)). 

The function M(a, b; Q is given by a convergent Taylor series 

Mia K) = 1 + !( + ^ + !MM^ + ... 
[ ' ' U + 61! + 6(6 + 1) 2! + 6(6+l)(6 + 2) 3! + 

As it is readily seen from the above series, M(a,b;Q is analytic in (a, 6,0 € C x {6 : \b — 1| < 
1/2} x C. For -a £ Z+, C > 0, 6 / 1, the function C/(a, 6; C) is defined by 

^/ \ t t / , *n 7r / r(a) , .i i,M(l + a — 6, 2 — 6; ON ^ 

r( " )f/ <°' fc « " S« ( r(i + a -V(t) M "- " : " c r ( 2- t) ) ' (5 ' 3) 

We assume that arg C = 0; this fixes the branch of 0~ 6 - The function T(a)U (a, 6; is meromorphic 
in a G C with poles at a = 0, —1, —2, . . . which correspond to the Landau levels — see (5.1). 

The r.h.s. of (5.3) is analytic in 6 with a removable singularity at 6 = 1; the limit as 6 — > 1 is 
easy to compute: 

T(a)U(a, 1; C) = -2Af£(a, 1; C) - M>, 1; C) - (2 7 + ^(o) - log C)M(a, 1; C), (5.4) 

where = M' b = Off, ^(a) = T'(a)/T(a) is the digamma function, 7 is Euler's constant 
7 = -V>(1) « 0.577, and log ( G M, ( > 0. 

Using the reflection formula for the ip function, 

tp(a) = ip(l — a) — 7rcot(7ra), 

let us rearrange formula (5.4) as 

T(a)U(a, 1; C) = M(a, Q + vr cot(vra)M(a, 1; Q, (5.5) 
M(a, C) = -2M' b {a, 1; C) - M>, 1; C) - (2 7 + Y>(1 - a) + log ()M(a, 1; C). (5.6) 

The function M(a, C) is analytic in a at the points —a G Z+. The singularities of r(a)[/(a, 1; C) 
written in the form (5.5) are easy to analyze, as they are due to the elementary function cot7ra. 
Incidentally, (5.5) gives a formula for the residues of the resolvent R(z); due to the identity 
M(—n, 1;0 = L n ((), this formula agrees with (3.1). 

Proof of Lemma 4-2. Substituting formula (5.1) into the integrals (4.3) and using (5.5), we see 
that in order to obtain the required asymptotic expansions, we need to analyse the asymptotics 
of the integrals 

/ (a + n) k (cot 7r a) u G (a) da, n -> 00, k G Z+, m£N, (5.7) 
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where 7„ is a positively oriented circle around a = —n with the radius 1/2, and G{a) is the analytic 
function 

G(a) = / F(xi, . . . ,Xj)~nl =1 M p (a, %\x p+ \ - x p | 2 )dxi • • • dxj, Xj + \ = x\. (5.8) 

Here 

/ i ^ \ 

F(xi . ..Xj) = V(xi) ■ ■ ■ V(xj) exp -f ^|x p+ i - x p \ 2 + if ^[x p ,x p+ i] , x j+ i = xi, 

V p =i p =i / 

each of the functions M p (a, () is either M(a, 1;£) or M(a, £) and at least one of the functions 
M p {a,() isM(a,l;C). 

Applying the residue formula to the integral (5.7), we see that the required statement follows 
from Lemma 5.2 below. ■ 

Lemma 5.2. For any F £ Cg°(R •?), Zet G(a) 6e given 6t/ (5.8), where each of the functions 
M p (a, Q is either M (a, 1; £) or M(o,() and ai Zeasi one of the functions M. p (a, () is M(a, 1;£)- 
TTten i/ie function G(a) and all of its derivatives G^ s \a), s ^ 1, admit asymptotic expansions in 
integer powers of | <x| x / 2 as a — > — oo ; a£i. 

Proof. The proof is based on using suitable asymptotic expansions of the functions M(a, 1; £) and 
M(a, in terms of Bessel functions and on application of Lemma 6.1. 

(i) First consider the special case when Ai p (a, Q = M(a, 1; Q for all p in (5.8). Our main tool 
will be a convergent expansion of M(a, b; () in terms of Bessel functions due to Tricomi [13]. For 
our purposes it suffices to consider the following range of parameters: 

Rea^-1, |Imo| < 1, |6-l|<l/2, < C < R (5.9) 

for some fixed R > 0. The expansion of [13] (see also [1, (13.3.7)]) reads: 

1 2 ) E(^T^) An^-i + m(V(26-4a)C), (5.10) 

' m=0 ^ ' 

where Jfe_i +m are Bessel functions and A m = A m (a, b) are the coefficients in the Taylor expansion 



m=0 

Note that Re (6 — 2a) ^ 1 and the principal values of ^ ( b ~^ a K ^ ^ ^ and (26 — 4a) are taken 

in (5.10). Due to a fast decay of the Bessel function J v (z) for v — > 00, the series (5.10) converges 
absolutely for the range of parameters (5.9). Take 6=1 and for a given N G N, write Tricomi's 
expansion (5.10) as 

M(o,l;C) =M^(a,() + M^(a,C), (5.12) 

where 

N-l , , xm/2 

i< 0) (a, C) = e^ 2 J] A m (a, 1) ^— J m (V(2 - 4a)C). 

m=0 ^ ' 
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Let us recall the argument of [13] which gives the estimate for M^\a,(). By inspecting the 
integral representation for the Bessel function, one obtains a uniform estimate 



|J m (V(2 - 4o)C)| < C for m£Z + , Rea<-1, |Ima| < 1, < C < R. (5.13) 

Next, one needs to estimate for the coefficients A m of the expansion (5.10). Applying Cauchy's 
theorem to the Taylor expansion (5.11) yields 



\A m \ ^r- m max|/(z) 

\z\=r 



for any r € (0, 1). Note that 



Choosing r = \a\~^ (one can take r = |a|~ 5 for any 1/3 < 5 < 1/2), one obtains for all sufficiently 
large \a\: 

\A m \ < C\a\^ m (l + C\a\-£y Rca ^ 2C\a\^ m . (5.14) 
Combining (5.13) and (5.14) gives for all sufficiently large \a\: 

oo r> rn/2 

\M^\a,0\^CY,W~' 1 



2 -4a 



0(\a\~T2 N ), Rea^-oo, |Ima| < 1. (5.15) 



m=N 

In the same way, the estimates (5.13) and (5.14) show that 

\M^\a,Q\ = 0(1), Rea^-oo, |Ima| < 1. (5.16) 
Substituting (5.12) into (5.8), we obtain 

G(a) = G^(a) + G^(a), 

where 

G w( a ) = / F(x 1 ,...,Xj)U 3 v=1 M^(a,§\x p+1 -x p \ 2 )dx 1 ---dxj, x j+1 = x\. 
By (5.16), (5.15), we get 

G$(a) = 0(\a\~& N ), Reo^-oo, |Ima| < 1. 
By Cauchy's formula for the derivatives, this entails 



d 



-T- ] G$(a) = 0(\a\-v N ), a^-oo, a G R. 
da J 

As N can be taken arbitrary large, we see that it suffices to prove that for any N > 0, all 
derivatives (^) S G$(a), s £ Z+, have an asymptotic expansion for a — > — oo, a£l. 

From (5.11) it follows that the coefficients A m (a, b) are polynomials in a and b. This observation 
reduces the problem to justifying the asymptotic expansion of the integral (5.8), where each of 
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the functions M p (a, () is ( m/2 J m (^(2 - 4a)() with some m G Z + . Such an expansion is provided 
by Lemma 6.1. 

(ii) Consider the general case. First let us obtain an expansion for M(a, Q similar to (5.10). 
Substituting (5.10) into the r.h.s. of (5.6), after a rearrangement we obtain 

i 00 / /■ \m/2 

M(a, C) = -(V(l - a) - log(- - a))e</ 2 £ A m ^— J m (V(2 - 4a)C) 

m=0 ^ ' 

- 2e^ 2 A ™ jm(V(2-4a)C) - e C/2 £ 5 ™ ^ J m (V(2-4a)C), 

m=0 ^ ' m=0 ^ ' 

(5.17) 



where 



do oa J av 



A fast decay of J m and J m as m — ► oo ensures convergence of the series and validates differentiation 
with respect to a and b. 

Using the expansion (5.17), we can complete the argument by following the same steps as 
in part (i) of the proof. First we need to obtain estimates for the remainders of the series in 
the r.h.s. of (5.17) in the strip Re a ^ —1, |Ima| ^ 1 (cf. (5.15)). The estimate for the 
remainder term of the first series in the r.h.s. of (5.17) is provided by (5.15). The estimate for 
the second series in the r.h.s. of (5.17) is obtained in exactly the same way by using the estimates 
| j (V(2 " 4a)C)| < C + C|log|(2 - 4o)C|| and (see [1, (9.1.22)]) 



\J m (V( 2 - 4a )C)l «S C, men, Rea^-1, |Ima| ^ 1, < ( ^ R 

instead of (5.13). In order to estimate the remainder term of the third series, we need an estimate 
on the coefficients B m . The coefficients B m are readily seen to be Taylor coefficients of the function 
(cf. (5.11)) 

g(z) = ^(z) + |£(z)) = f(z) log(l - z 2 )- 1 / 2 
which similarly to (5.14) gives 

\B m \ < C\a\T2 m . 

This gives the analogue of the estimate (5.15) for the third series in the r.h.s. of (5.17). 

Next, the function if)(l — a) — log(| — a) in (5.17) admits asymptotic expansion in integer 
powers of a -1 as Re a — > — oo (see [1, 6.3.18]). Thus, we have reduced the problem to jus- 
tifying an asymptotic expansion of the integral (5.8), where each of the functions Ai p (a, £) is 
either ( m / 2 J m (y/(2 - 4a) Q or ( m / 2 j m (y/(2 - 4a) Q and at least one of the functions M p (a, Q is 
Q m / 2 J m {y/{t^4ajC,). Finally, the formula [1, (9.1.66)] 



m— 1 



(C/2)* 



reduces the problem to Lemma 6.1. 
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6 Asymptotics of integrals containing Bessel functions 

Lemma 6.1. Let F £ (R 2j ) . Define a function Gi(n), k > 0, by 

Gi(k) := / F(x 1 , . . .Xj)T\ J m (K\x p - x p+ i\)\x p - x p+ i\ mp dxi- ■ -dxj, x j+1 = Xj (6.1) 

where each of the functions J mp (£) is either J mp (C) or ^m p (C) with some m p £ Z+, and at least 
one of the functions J mp (() is J mp {()- Then the function G\{k) and all of its derivatives g[ s \k), 
s ^ 1, have a complete asymptotic expansion in integer powers of n~ l for n — ► +oo. 

Proof. Recall the identity [1, (9.1.27)] 

r +1 J,+i(C) = 2"CMQ - C 2 (C"" 1 ^-i(C)), !«/ = Ju or J, = Y v . 

This identity allows us to reduce the problem to the case when all the indices m p in the integral 
(6.1) are or 1. Next, assume for the convenience of notation that the first / functions J in the 
integral (6.1) are the Neumann functions Y, and the remaining j — I functions are the Bessel 
functions J. For n\ > 0, k 2 > 0, . . . , kj > define 

G 2 (Kl,...,Kj) 



f 



F(x±,...Xj) Y\_Y (k p \x p - x p+1 \) 'Yl Jo(k p \x p - x p+ i\)dxi ■ ■ -dxj, x j+1 = x 1 . (6.2) 
p=i p=i+i 



Formulae 

^M = -|x|j 1 («N) > d ^^ = -\x\YMA) 

show that it suffices to obtain an asymptotic expansion of the functions 

{llY 1 ■ ■ ■ (4)* ■ ■ ■ «*) Ui A, e z + (6.3) 

for k — ► oo. We shall obtain an asymptotic expansion for the function Gs(k) := G 2 {k, ■ ■ ■ , k) as 
k — > oo. From the construction it will be clear that the derivatives (6.3) can be dealt with in the 
same way. Let us make a change of variables in the integral (6.2). Denote 

y p = x p -x p+1 , p=\,...,j-l, z = xi H \-Xj, 

F 1 {y) = [ F( Xl . . . x 3 )dz, F 1 G C^°(E 2j, ~ 2 ). 

Then we obtain 

, i j'-i 

G 3 (k)= / Fi(i/)Jo(K|j/i + --- + y j _i|)TT y o(«|y P |) IT J o(«|l/ P |)di/. (6.4) 

Next, we use the formulae 

If If piK-uy 

J (K\y\) = - e iKu yd(u 2 -l)du, Y (K\y\) = - -v.p. / - d«. 
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Substituting these formulae into (6.4), we obtain 



G 3 (k) = v.p. / 



8{v? l+1 -l)---5{u)-l) 
(u? - 1) • ■ ■ (u? - 1) 



F 2 (k(Ui — Uj,U2 — Uj, . . . , Uj-l — Uj))du\ . . . dUj, 



where F 2 is (up to a multiplicative constant) the Fourier transform of F\. Note that F2 belongs 
to the Schwartz class S(R 2j ~ 2 ). 

In order to simplify the last integral, we introduce some notation. Let us use the polar 
coordinates U{ = r l J 2 u)i, where tDj = (coscjj, sincjj) G M 2 , oj-i G T = M/27rZ. Denote also u> = 
(ui, . . . ,u>j-i) G TP -1 , r = (n, . . . ,77) G R+. Define the function 



f(r,u,Uj) = {r\ /2 ui - Uj, . . .,rj /2 Ui - - tij, . . - toj) G 

With this notation, we have 

F 2 (nf(r,u,Uj)) 



j-2 



Gs(k) = 2 J / dcjj / do; v.p. / dr 

Jt Jv- 1 Jm!, 



(n - 1) ■ ■ ■ (r, - 1) ' 



(6.5) 



Note that 



f(r,u,Uj) =0& (r = (!,...,!) €R l + and u= {uij, ...,uj)e TV -1 ) 



to j, ... ,ujj). Let us show 



and rank/'(r, w, Wj) = Z + j — 1 at the point r = (1,...,1), uj 
that only an arbitrary small neighbourhood of the point r = (1, . . . , 1), uj = (uij, . . . ,u>j) gives 
contribution to the asymptotics of the integral (6.5). First recall some estimates for the principal 
value integrals. Let 5 > and (j) G C°°(-8, 5). Then 



v.p. 



-5 X 

and using the Cauchy-Schwartz inequality, one obtains 



f 5 , 

dx = <fi (x) log\S/x\dx, 
J-s 



v.p. J ' dx 



-5 x 



C{5)\W\\ L > { -5,8y 



(6.6) 



Similarly, for d> G C°° ([-<$, 5} 1 ) 



v.p. 



-dx 



(-5,5)1 Xi---Xi 



dx\ ■ ■ ■ dxi 



LH(-8,8Y) 



(6.7) 



Denote U = (1 — e, 1 + e) 1 x (uij - e, ujj + e) j 1 C R l + x TP 1 where e > is sufficiently small. Let 
us show that 

G 3 (k) = 2"' / c%v.p. / dr^ -^^p^p^ +0(k-°°). (6.8) 
7t 7c/ (ri - 1) ■ ■ ■ (n - 1) 

For simplicity consider the case j = 2, I = 1. Then for k — > 00 one has 

dr\ 



f duj 2 f du>i f dri F 2 (Kf(r 1 ,u 1 ,u) 2 )) 

Jt Jt JR+\(i-e,i+e) n - 1 



sup |F 2 (K/(n,u;i,W2))| = 0(«-°°), 

|ri-l|>e 
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and also by (6.6) 



duj 2 v.p. 



l+e 



dri 



l-e H — 1 Jt\(uj 2 -s,u> 2 +s) 



duj 1 F 2 (Kf(r 1 ,u 1 ,uj 2 )) 

f 

Ji- 



1+e 

duii I dr\ 



dF 2 {Kf{r 1 ,0Ji,uj2)) 



dr± 



= 0(k 



as F 2 is the Schwartz class function. 

Thus, it suffices to prove an asymptotic expansion of the integral in the r.h.s. of (6.8). The 
asymptotic expansion of the integral over (r, cj) is provided by Lemma 6.2 below. It remains to 
prove that the asymptotic expansion given by Lemma 6.2 is uniform in uij. In order to show this, 
let us introduce the matrix 

'cos r — sin r\ 

, t G T, 



sin r cos r 



and the vector e = (1, 0) G M 2 . We rewrite the function f(r,uj,ujj) in the form 

f(r, u, ujj) = ((r{ /2 VVi - !)e, • • • , {r] ,2 il> Vl - l)e, - l)e, . . . , (^_i - l)e) 

= ^wjf(r,i/,0), u p = u p -u; j , v = (v 1 , . . . , G T- 7-1 , (6.9) 

where ^ T = diag(Vv, • • • , Vv) is a matrix in IR 2 - 7 " 2 . Substituting (6.9) into (6.8), we obtain 

-if [ F 2 (K^/(r,z/,0)) 

G 3 {k) = 2 J / dujjv.p. / drdz;- — ^ — — 

Jt Jc/i (ri - 1) • • • (r, - 1) 



, t/i = (l-e, l+e)' x(-e,< 



C Mi xT 



(6.10) 



From the last formula and the proof of Lemma 6.2 it is clear that the expansion given by Lemma 6.2 
is uniform in coj; integrating this expansion over wj, we get the required expansion for Gs(k). ■ 

Lemma 6.2. Let f G C°°(R m , R n ), m^n, and suppose that /(0) = 0, rank /'(()) = m. Then for 
any sufficiently small open neighbourhood of zero U C M m , any F G 5(M n ), and I G {0, 1, . . . , m}, 
the integral 

f F{Kf{x)) , 
IU) = v.p. / v Jy " dx 
Ju xi---xi 

has a complete asymptotic expansion 

/(«) = K l ~ m (co + ci^T 1 + c 2 kT 2 H ), K — ► oo. 

Proof. Choose U sufficiently small so that 

\f(x)-f'(0)x\^±\f'(0)x\, VxeU. 
For a given N G N, N > /, let us prove that 

N-m-l 

J-m „ i r^f.J-N\ 



I{K) = K l ~ m ^ + °( k1 ~ 

i=0 



K — > OO. 



(6.11) 



(6.12) 
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N 



By Taylor's formula for f(x) and F{nf(x)), 

f(x) = f'(0)x + f 2 (x), f 2 (x) = £ -,f is) (0)x s + f N (x), (6.13) 

8 = 1 S ' 

N 

F( K f(x)) = F( K f\0)x + K f 2 (x))=J2 F{q \^f\^)^f2(x)) q + F N (x,K), (6.14) 

q=0 

Fn(x, K ) = ± j\l- t ) n (£) N+1 F(Kf'(0)x + Krf 2 (x))dr. (6.15) 

Here we use simplified notation; f( s \0)x s stands for the polylinear form of the s'th differential of 
/ at zero, etc. 

Substituting (6.13) into (6.14) and collecting the terms that contain and that do not contain 
/n(x) into two different sums, we can write 



F( K f(x))= F q (Kx)P q (x)+ Y F q (Kx)g q (x)+F N (x,K). 



(6.16) 



Here both sums over q are finite, F q € S(M m ) are obtained from various components of derivatives 
of F, P q (x) are polynomials in x of degree q, and g q € C°°(U) are functions, satisfying the 
polynomial type estimates 



a 

dx, 



01 



9q(x) 



/?! + ••• + 0i < q, xeU. 



(6.17) 



Consider the terms obtained by substitution of the r.h.s. of (6.16) into the integral (6.10). First, 
using the estimate (6.7) and the fact that F q is a Schwartz class function, we obtain 



f F q (nx) f 
.p. / — F q {x)dx = v.p. / 

Ju x i ' ' ■ x i Jw 



Fq{Kx) P q (x)dx + 0(k-°°) 
x\---xi 



K l- q -m wv 



J 



Fq(x) 
Xl-'-Xl 



P q {x)dx + 0{n °°), K — > oo. 



So, these terms will give contribution to the asymptotics (6.12). 

Next, consider the terms obtained by substitution of the second sum in (6.16) into the integral 
(6.10). Using (6.17), we obtain the estimate 



d\F q {Kx)g q {x)) 



dx\ ■ ■ ■ dxi 



L 2 (U) 



By (6.7), it follows that all the corresponding integrals are 0(k 1 n ) as k — > oo. 
Finally, consider the term F^(x, k). By (6.11), we obtain for some c > 0: 

\f'(0)x + Tf 2 (x)\ ^ ^\f'(0)x\ >c\x\, xeU, re (0,1). 



Using this fact, we obtain 



\F n (x,k)\^Ck n+1 sup \F( N+1 \y)f 2 (x) N+1 \ 

\y\^nc\x\ 
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and therefore 

\\FN(;K)\\mu) < 0{k-^- n - v ), k - +00. 
Similarly, one can prove the estimate 

8 1 Fn(x, k) 



dx\ ■ ■ ■ dxi 



L2(C/) 



By (6.7), it follows that the integral of Fn is 0(k 1 n 1 ) and will only give contribution to the 
remainder term in (6.12). ■ 
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